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Abstract 

The bulk of the hydrogen in the universe transformed 
from neutral to ionized somewhere in the redshift interval 
5 ^ ^ ^ 40, most likely due to ionizing photons produced 
by an early generation of stars or mini-quasars. The re- 
sulting free electrons, interacting with the CMB photons 
via Thomson-scattering, are a mixed blessing, providing 
both a probe of the epoch of the first stars and a contam- 
inant to the pristine primary anisotropy. Here we review 
our current knowledge of reionization with emphasis on 
inhomogeneities and describe the possible connections to 
CMB anisotropy. 

Subject headings: Reionization, cosmic microwave background, 
foregrounds 

1. Introduction 

One of the most remarkable observational results in 
the last three decades in cosmology is the lack of the 
Gunn-Peterson trough in the spectra of high-redshift 
quasars and galaxies. This finding implies that the in- 
tergalactic medium (IGM) is highly ionized at least out 
to z ~ 5, the redshift of the most distant known quasars 
and galaxies]^. Since the primeval plasma recombined at 
z ^ 1100, its subsequent ionization requires some form of 
energy injection into the IGM, naturally attributable to 
astrophysical sources. The two most popular examples 
of such sources are an early generation of stars, residing 
in sub-galactic size clusters (hereafter "mini-galaxies"), 
or accreting massive (~ lO^M©) black holes in small 
( ^ IO^Mq) halos (hereafter "mini-quasars"). 

Reionization is interesting for a variety of reasons. 
First, the most fundamental questions are still unan- 
swered; What type of sources caused reionization, and 
around what redshift did it occur? How were the reion- 
izing sources distributed relative to the gas? What 



^Currently the best lower limit on the reionization redshift comes 
from the detection of high-z Ly« emission lines, implying z > 5.64 
(see discussion below). 



was the size, geometry, and topology of the ionized 
zones, and how did these properties evolve? Second, 
and more directly relevant to the topic of this review, 
reionization leaves distinctive signatures on the cosmic 
microwave background (CMB) through the interaction 
between the CMB photons and free electrons. Although 
it may complicate the extraction of cosmological param- 
eters from CMB data, this "contamination" could yield 
important information on the ionization history of the 
IGM. Through reionization, the CMB is a useful probe 
of nonlinear processes in the high-redshift universe. 

Thomson scattering from free electrons affects the 
CMB in several ways. Because the scattering leads 
to a blending of photons from initially different lines 
of sight, there is a damping of the primary tempera- 
ture anisotropy. On the other hand, a new secondary 
anisotropy is generated by what can be thought of as a 
Doppler effect: as photons scatter off free electrons, they 
pick up some of their peculiar momentum. Finally, the 
polarization dependence of the Thomson cross section 
creates new polarization from the initially anisotropic 
photon field. 

The inhomogeneity of reionization affects these pro- 
cesses in two different ways. First, the inhomogeneity 
of the medium affects how the mean ionization fraction 
evolves with time. This is important since the damping, 
Doppler and polarization effects all occur even in the ide- 
alized case of spatially homogeneous reionization. Sec- 
ondly, the spatial fluctuations in the ionization fraction 
can greatly enhance the contribution from the Doppler 
effect at small angular scales. Polarization and damping 
are much less influenced by the inhomogeneity, as we will 
see below. 

The aim of the present review is to summarize our 
current knowledge of reionization, describe the possible 
connections to the CMB, and assess what we can hope 
to learn about the high-redshift IGM from forthcoming 
observations of the CMB anisotropics. 

2. Reionization 

2.1. Homogeneous Reionization and the Reion- 
ization Redshift 

If the IGM was neutral, its optical depth to Lya ab- 
sorption would be exceedingly high: 
Tigm -- 105(r2b/i/0.03)[(l + z)/6]3/2, wiping out the flux 
from any source at observed wavelengths shorter than 
ALya(l + z). The spectra of high-redshift quasars and 
galaxies reveal Lya absorption by numerous discrete Lya 
forest clouds, separated in redshift, rather than the con- 
tinuous Gunn-Peterson (GP) trough expected from a 



neutral IGM. The detection of the continuum flux in- 
between the Lya clouds implies that Tigm is at most of or- 
der unity, or, equivalently, an upper limit on the average 
neutral fraction of a; = nni/^-tot ^ 10~^. The intergalac- 
tic medium is therefore highly ionized. Note that there 
is no sharp physical distinction between the intergalactic 
medium and the low column density (iVni ^ 10^^ cm~^) 
Lya forest (Reisenegger & Miralda-Escude 1995). These 
Lya clouds fill most of the volume of the universe, are 
highly ionized, and account for most of the baryons 
(implied by Big Bang nucleosynthesis) at high redshift 
(Weinberg et al. 1997; Ranch et al. 1997). 

The energy requirement of reionization can easily 
be satisfied either by the first mini-galaxies or mini- 
quasarsQ. Nuclear burning inside stars releases several 
MeV per hydrogen atom, and thin-disk accretion onto 
a Schwarzschild black hole releases ten times more en- 
ergy, while the ionization of a hydrogen atom requires 
only 13.6 eV. It is therefore sufficient to convert a frac- 
tion of < lO^'^ of the baryonic mass into either stars 
or black holes in order to ionize the rest of the baryons 
in the universe. In a homogeneous universe, the ratio 
of recombination time to Hubble time in the IGM is 
trec/tuuh ~ [(1 + z)/ll]~^^''^ SO that cach H atom would 
only recombine a few times at most until the present 
time. Unless reionization occurred at z ^ 10, recombi- 
nations therefore would not change the above conclusion, 
and the number of ionizing photons required per atom 
would be of order unity. Based on three-dimensional 
simulations, Gnedin & Ostriker (1997) argued that the 
dumpiness of the gas increases the average global recom- 
bination rate by a factor C =< p"^ > / < p >^^ 30 — 40 
(the averages are taken over the simulation box). This 
would imply a corresponding increase in the necessary 
number of ionizing photons per atom. However, in a more 
detailed picture of reionization in an inhomogeneous uni- 
verse (Gnedin 1998, Miralda-Escude et al. 1999), the 
high-density regions are ionized at a significantly later 
time than the low density regions. If this picture is cor- 
rect, then one ionizing photon per atom would suffice 
to reionize most of the volume (filled by the low-density 
gas) ; recombinations from the high-density regions would 
only contribute to the average recombination rate at a 
later redshift. 

The reionization redshift can be estimated from first 
principles. In the simplest picture [first discussed by 
Arons & Wingert (1972) in the context of quasars], a 
homogeneous neutral IGM is percolated by the Individ- 
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Fig. 1. — Schematic stages of reionization by early 
mini-galaxies. The small dots denote small halos with 
Tvir < lO^K, the larger dots correspond to larger halos 



with Tvir > lO'^K. Top panel 



40: The first ha- 



los form, both small and large halos create HII zones. 
Middle panel - z ~ 20: The UV background suppresses 
star- formation in small halos, only large halos can pro- 
duce HII zones. Bottom panel - z ~ 15: the ionized 
zones from the large halos overlap. 



^More exotic possibilities that we do not discuss here include pri- 
mordial black holes (e.g. Gibilisco 1996); cosmic rays (e.g. Nath 
& Bierman 1993); winds from supernovae (e.g. Ostriker & Cowie 
1981); or decaying neutrinos (Sciama 1993). 



ual HII regions growing around each isolated UV source. 
The universe is reionized when the ionized bubbles over- 
lap, i.e. when the filling factor of HII regions reaches 
unity^. Several authors studied this problem (Carr et 
al. 1984, Couchman & Rees 1986, Fukugita & Kawasaki 
1994, Meiksin & Madau 1993; Shapiro et al. 1994, 
Tegmark et al. 1994; Aghanim et al. 1996; Haiman & 
Loeb 1997). Specifically, Haiman & Loeb (1998a) have 
used the Press-Schechter formalism to describe the for- 
mation of halos that potentially host ionizing sources, 
in order to estimate the redshift at which overlap oc- 
curs. In this study, star and quasar black hole formation 
was allowed only inside halos that can cool efRciently 
via atomic hydrogen, since molecular hydrogen would be 
photodissociated earlier, as argued by Haiman, Abel & 
Rees (1999) and Haiman, Rees & Loeb (1997). This re- 
quirement translates to a minimum virial temperature of 
- lO^K, or halo masses above - 10®Mo[(l + z)/!!]-^/^ 
(see Fig. |l|) . For both types of sources, the total amount 
of light they produce was calibrated using data from red- 
shifts z <5. The efficiency of early star formation can be 
estimated from the observed metallicity of the intergalac- 
tic medium (Songaila & Cowie 1996, Tytler et al. 1995), 
utilizing the fact that these metals and the reionizing 
photons likely originate from the same stellar popula- 
tion. Similarly, the efficiency of black hole formation 
inside early mini-quasars can be constrained to match 
the subsequent evolution of the quasar luminosity func- 
tion at redshifts z < 5 (Pei 1995). These assumptions 
generically lead to a reionization redshift 8 ^ 2 ^ 15. 
The uncertainty in the redshift reflects a range of cos- 
mological parameters (i/o, i^m, i^b), CDM power spectra 
(n, as), and efficiency parameters for the production and 
escape fraction of ionizing photons. In the stellar reion- 
ization scenario, the predicted redshift could be above or 
below this range, if the initial mass function was strongly 
biased (relative to Scalo 1986) towards massive or low- 
mass stars. 

2.2. Importance and Treatments of Inhomogeneities 

Although the naive picture of overlapping HII spheres 
may give a good estimate for the "overlap redshift" , it is 
likely to be a crude and incomplete description of the 
real process of reionization. The ionizing sources are 
expected to be located in the highly overdense regions 
that form the complex large scale structure of the uni- 
verse - along filaments or sheets, or at their intersec- 
tions (see Fig. |2|). The ionizing radiation will have to 
cross the local dense structures before propagating into 




^The expansion of spherical ionization fronts around steady sources 
can be calculated analytically (Shapiro & Giroux 1987). 



Fig. 2. — A two-dimensional slice of log((5b), the bary- 
onic overdensity at redshift z=5 in a three-dimensional 
simulation (see Zhang et al. 1998). 



the IGM, which itself has significant density fluctuations. 
Several important consequences of these inhomogeneities 
can be noted. First, reionization could occur "outside-in" 
(Miralda-Escude et al. 1999), if the radiation from a typi- 
cal source escapes through a relatively narrow solid angle 
from its host, without causing significant local ionization. 
In this case, the low-density regions could be ionized 
significantly earlier than the high-density regions. An 
important effect in this case could be the shadowing of 
ionizing radiation by the dense concentrations within the 
IGM. The same type of shadowing effect was found to be 
important at low redshift. Based on the number of Lya 
absorbers, Madau, Haardt & Rees (1998) have concluded 
that as much as ~ 50% of the UV flux is absorbed by 
these systems. The redshift at which the IGM becomes 
optically thin to the ionizing continuum is therefore sig- 
nificantly delayed relative to the overlap epoch, when 
the ionized zones first percolate. An alternative possibil- 
ity is that the ionizing sources are densely spaced along 
filaments, and their radiation ionizes most of the dense 
filament before escaping into the IGM. In this "inside- 
out" case, the description of ionized zones surrounding 
their host source would be more applicable. An impor- 
tant feedback during and after reionization is that the 
newly established UV background can photo-evaporate 
the gas from small (wcirc ^ 10) halos (Shapiro, Raga & 
Mellema 1997; Barkana & Loeb 1999). A related point, 
important for the detection of the GP trough, is that 
at the time of overlap of the HII zones, each individual 
HII region could still have a non-negligible Lya optical 



depth, due to its residual HI fraction. The GP trough 
therefore disappears only from the spectra of individual 
sources located below the later redshift when the average 
Lya optical depth of the IGM drops below unity (Haiman 
& Loeb 1998b; Miralda-Escude et al. 1999; Shapiro et al. 
1987). Finally, the inhomogeneity of the reionized IGM 
leads to generation of anisotropics in the CMB which are 
strongly suppressed in the homogeneous case, as we shall 
see below. 

Reionization is now also being addressed by three- 
dimensional simulations, which can eventually shed light 
on some of the above issues. The first such studies ap- 
proximated radiative transfer by assuming an isotropic 
radiation field, but were able to take into account the in- 
homogeneities in the gas distribution (Ostriker & Gnedin 
1996; Gnedin & Ostriker 1997). The resulting reioniza- 
tion redshift of z ~ 7 in a ACDM model is in good agree- 
ment with the semi-analytical estimates, given the dif- 
ferences in the assumed cosmology, power spectra, and 
star formation efficiency. However, an isotropic radia- 
tion field is necessarily a crude approximation, at least 
in the beginning stages of reionization, when each source 
still has its own isolated ionized zone. In order to un- 
derstand the progressive overlap of ionized zones, includ- 
ing effects such as self-shielding and shadowing by dense 
clumps, it is necessary to incorporate radiative transfer 
into three dimensional simulations (Norman et al. 1998). 
The numerical algorithms for the first such attempts are 
described in Abel, Norman & Madau (1998) and Ra- 
zoumov & Scott (1998). Results from an approximate 
treatment (Gnedin 1998) also support the outside~in pic- 
ture of reionization described above. 

2.3. The Nature of Reionizing Sources 

Several authors have emphasized that the known pop- 
ulation of quasars or galaxies provide ^10 times fewer 
ionizing photons than necessary to reionize the universe 
(see, e.g. Shapiro et al. 1994, or Madau et al. 1999, and 
references therein). The two most natural candidates 
for undetected ionizing sources are the mini-galaxies 
and mini-quasars expected to be associated with small 
(~ 10^~^° M0) dark halos at 2: ~ 10. As there are no 
compelling ab initio theoretical arguments to strongly fa- 
vor one type of source over the other, the best hope is to 
distinguish these two possibilities from observations. One 
expects at least three major differences between mini- 
galaxies and mini-quasars: their spectra, their absolute 
brightness (or, alternatively, space density), and the an- 
gular size of their luminous regions. 

The fiux from stellar populations drops rapidly with 
frequency above the ionization threshold of hydrogen, 
while the spectra of mini-quasars are expected to be 



harder and extend into the X-ray regime. Stars could 
therefore not reionize Hell, while quasars with typical 
spectra (Elvis et al. 1994) would reionize HcII at approx- 
imately the same redshift as H (Haiman & Loeb 1998a). 
However, recent high-resolution spectra of z ~ 3 quasars 
have shown a widely fluctuating Hell/HI optical depth 
(Reimers et al. 1997). The authors interpreted this ob- 
servation as evidence for HcIII regions embedded in an 
otherwise Hell medium, i.e. a detection of the Hell reion- 
ization epoch (see also Wadsley et al. 1998). If this in- 
terpretation is verified by future data, it would constrain 
the number of mini-quasars with hard spectra extend- 
ing to X-rays at high redshifts (Haiman & Loeb 1998b; 
Miralda-Escude 1998; Miralda-Escude & Rees 1994). At 
present, a plausible alternative interpretation of the Hell 
observations is that the observed Hell optical depth fiuc- 
tuations are caused by statistical fluctuations in the IGM 
density or the ionizing background flux (Miralda-Escude 
et al. 1999), rather than by the patchy structure of 
Hell/HelH zones. 

The X-ray background (XRB; Miyaji et al. 1998; 
Fabian & Barcons 1992) might provide another use- 
ful constraint on mini-quasar models with hard spectra 
(Haiman & Loeb 1998c). These models ovcrpredict the 
unresolved flux by a factor of ^ 2 — 7 in the 0.1-1 keV 
range. If an even larger fraction of the XRB will be re- 
solved into low-redshift AGNs in the future, then the 
XRB could be used to place more stringent constraints 
on the X-ray spectrum or the abundance of the mini- 
quasars. 

A second distinction between galaxies and quasars can 
be made based on their absolute luminosity and space 
density. On average, based on their z < 5 counterparts, 
high-redshift mini-quasars are expected to be roughly 
~ 100 times brighter, but ~ 100 times more rare, than 
galaxies. At the overlap epoch, HII regions from mini- 
quasars would therefore be larger but fewer than those 
from mini-galaxies. The typical size and abundance of 
HII regions has at least three important implications. 
First, it affects the CMB anisotropics; the larger the HII 
regions, the larger the effect on the CMB (Gruzinov and 
Hu 1998, Knox et al. 1998), as we will see below. Sec- 
ond, large enough HII regions ( ^ IMpc) would allow 
gaps to open up in the GP trough; the flux transmit- 
ted through such gaps around very bright quasars could 
be detectable (Miralda-Escude et al. 1999). Third, if 
sources form inside the rare, high-a peaks, they would 
be strongly clustered (Knox et al. 1998). This clustering 
would increase the effective luminosity of each source, 
and would enhance both of the effects above. 

Finally, a constraint on the number of mini-quasars 
at z > 3.5 can be derived from the Hubble Deep Field 



(HDF). Mini-quasars are expected to appear as faint 
point-sources in the HDF, unless their underlying ex- 
tended host galaxies are resolved. The properties of faint 
extended sources found in the Hubble Deep Field (HDF) 
agree with detailed semi-analytic models of galaxy for- 
mation (Baugh et al. 1998). On the other hand, the 
HDF has revealed only a handful of faint unresolved 
sources, and none with the colors expected for high red- 
shift quasars (Conti et al. 1999). The simplest mini- 
quasar models predict the existence of ~ 10 — 15 B- 
band "dropouts" in the HDF, inconsistent with the lack 
of detection of such dropouts up to the ~ 50% com- 
pleteness limit at F « 29. To reconcile the models 
with the data, a mechanism is needed that suppresses 
the formation of quasars in halos with circular veloci- 
ties Wcirc ^ 50 — 75 km s~^ (Haiman, Madau & Loeb 
1999). This suppression naturally arises due to the 
photo-ionization heating of the intergalactic gas by the 
UV background after reionization. Forthcoming data on 
point-sources from NICMOS observations of the HDF 
(see Thompson et al. 1998) could further improve these 
constraints. 

2.4. When vi^as the Universe Reionized? 

As described above, theoretical expectations place the 
reionization redshift at 8 ^ z ^ 15. What about observa- 
tions? At present, the best lower limit in the reionization 
redshift comes from the detection of high-redshift Lya 
emitters (Hu et al. 1998; Weymann et al. 1998). As 
argued by Miralda-Escude (1998), the damping wing of 
Lya absorption by a neutral IGM has a large residual 
optical depth that would severely damp any Lya emis- 
sion line. Currently the highest redshift at which a Lya 
emitter is seen is z = 5.64; the existence of this object im- 
plies that reionization occurred prior to this redshift (see 
Haiman & Spaans 1998). Note that HII regions around 
individual sources, excepting only the brightest quasars, 
would be too small to allow the escape of Lya photons, 
because the damping wings of the GP troughs around 
the HII region would still overlap. 

The best upper limit on the reionization redshift can 
be obtained from the CMB anisotropy data. Several ex- 
periments have revealed a rise in the power on small an- 
gular scales, and a drop at even smaller scales (e.g., Bond 
et al. 1999), showing evidence for the first Doppler peak 
expected from acoustic oscillations in the baryon-photon 
fluid prior to recombination (e.g., Hu & Sugiyama 1994). 
These observations can also be used to set limits on the 
electron scattering optical depth (or the corresponding 
reionization redshift) that would suppress the Doppler 
peak. From a compilation of all the existing measure- 
ments, Griffiths et al. (1998) have derived the stringent 



(although model-dependent) constraint z ^ 40 on the 
reionization redshift - if reionization occurred earlier, the 
electron scattering optical depth would have reduced the 
amplitude of the anisotropics below the observed level. 
Note that another upper limit results from the spectral 
distortion of the CMB caused by scattering on the reion- 
ized IGM, which is constrained by the upper limit on 
the Compton y-parameter measured by COBE. The in- 
ferred upper limit on the reionization redshift, however, 
is much weaker, z <J 400 for typical parameters in a low- 
density universe (Griffiths et al. 1998, see also Stebbins 
& Silk 1986). In summary, present observations have 
narrowed down the possible redshift of the reionization 
epoch to 6 ^ z ;$ 40, a relatively narrow range that is 
in good agreement with the theoretical predictions de- 
scribed above. 

2.5. Future Observational Signatures 

Further observational progress in probing the reioniza- 
tion epoch could come either from the Next Generation 
Space Telescope (NGST), or from more precise measure- 
ments of the CMB anisotropics by Myli^ and P/anc^. 
NGST, scheduled for launch in 2007, is expected to reach 
the ^ InJy sensitivity^ required to detect individual 
sources to z ~ 10, and to perform medium-resolution 
spectroscopy to z ^ 8. If reionization occurred close to 
the low end of the allowed redshift range, 6 ^ z ^ 8, then 
it may be possible to infer the redshift directly from the 
spectra of bright sources. One specific method relies on 
the spectrum of a bright source just beyond the reion- 
ization redshift, so that the individual Lya, Ly/?, and 
other GP troughs do not overlap in frequency, leaving 
gaps of transmitted flux (Haiman & Loeb 1998b). The 
measurement of this transmitted flux would be possible 
with NGST, despite absorption by the high-redshift Lya 
forest, as long as the number density of absorbers does 
not rise much more steeply with redshift than an extrap- 
olation of the current z < 5 data would imply (Fardal et 
al. 1998). If reionization is gradual, rather than abrupt, 
than this method would measure the redshift at which 
the GP optical depth drops to near unity, i.e. the final 
stages of the reionization epoch. An alternative signa- 
ture to look for would be the precise shape of the damp- 
ing wing of Lya absorption from the neutral IGM along 
the line of sight to the source (Miralda-Escude 1998). 
The shape, if measured by high-resolution spectroscopy, 
could be used to determine the total optical depth of the 
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IGM. The caveat of this method is the inabihty to dis- 
tinguish the neutral IGM from a damped Lya absorber 
along the line of sight, close to the source. 

An alternative signature is the background Lya emis- 
sion from the reionization epoch. Recombinations are 
slow both at high rcdshifts, when the IGM is still neu- 
tral, and at low rcdshifts when the IGM density is low. 
As shown by simulations (Gnedin & Ostriker 1997), the 
global recombination rate has a pronounced peak around 
the reionization epoch. The resulting recombinant Lya 
emission has been computed by Baltz, Gnedin & Silk 
(1998), and a detailed detectability study has shown that 
this signal could be measured by NGST, or perhaps even 
by HST (Shaver et al. 1999). 

Yet another signature could result from the 21 cm hy- 
perfine transitions in the IGM before reionization. Prior 
to reionization, the excitation of the 21 cm line in neu- 
tral HI atoms depends on the spin temperature. The 
coupling of the spin temperature to that of the CMB is 
determined by the local gas density, temperature, and 
the radiation background from the first mini-galaxies 
and mini-quasars. In general, the 21 cm line could be 
seen either in absorption or emission against the CMB, 
and could serve as a 'tomographic' tool to diagnose the 
density and temperature of the high redshift neutral gas 
(see, e.g. Madau et al. 1997; Scott & Rees 1990). If 
reionization occurred at 6 ^ 2; ^ 10, the redshifted 21 
cm signals would be detectable by the Giant Metrewave 
Radio Telescope; a study of the effect at higher rcdshifts 
would be possible with next generation instruments such 
as the THousand Element Array, or the Square Kilome- 
ter Array (Shaver et al. 1999). 

The reionization redshift could turn out to be closer to 
the high end of the allowed range, 10 ^ z ^ 40, rendering 
direct detection of emission from the reionization epoch 
implausible (except perhaps the 21cm signal). There is, 
however, a fortunate "complementarity", since the elec- 
tron scattering optical depth increases with reionization 
redshift. Due to the effect on the polarization power 
spectrum on large angular scales, MAP and Planck may 
be able to discern an electron scattering optical depth as 
small as a few percent (Zaldarriaga et al. 1997, Eisen- 
stein et al. 1998, Prunet et al. 1998). By temperature 
anisotropy data alone, they may be able to determine 
an optical depth as small as ^ 20% due to the damping 
effect mentioned earlier — and to be discussed in more de- 
tail later. 

The reionizing sources may also change the spectral 
shape of the CMB. The dust that is inevitably pro- 
duced by the first type II supernovae, absorbs the UV 
emission from early stars and mini-quasars and re-emits 
this energy at longer wavelengths. Loeb & Haiman 



(1997) have quantified the resulting spectral distortion 
in Press-Schechter type models, assuming that each type 
II supernova in a Scalo IMF yields O.SM© of dust with 
the wavelength-dependent opacity of Galactic dust, uni- 
formly distributed throughout the intergalactic medium. 
Under these assumptions, the dust remains cold (close 
to the CMB temperature), and its emission peaks near 
the CMB peak. The resulting spectral distortion can be 
expressed as a Compton y-parameter ~ 10^"'', near the 
upper limit derived from measurement of the CMB spec- 
trum by COBE (Fixsen et al. 1996). A substantial frac- 
tion (^ 10-50%) of this total y-parameter results sim- 
ply from the direct far-infrared emission by early mini- 
quasars and could be present even in the absence of any 
intergalactic dust. 

Inhomogeneities in the dust distribution could change 
these conclusions. Instead of being homogeneously mixed 
into the IGM, the high-redshift dust may remain concen- 
trated inside or around the galaxies where it is produced. 
In this case, the average dust particle would see a higher 
flux than assumed in Loeb & Haiman (1997), so that 
the dust temperature would be higher, and dust emis- 
sion would peak at a shorter wavelength. The magni- 
tude of the spectral distortion would be enhanced, and 
may account for the recently discovered cosmic infrared 
background (Pugct et al. 1996, Schlegel et al. 1998, 
Fixsen et al. 1998, Hauser et al. 1998) as shown by the 
semi-analytic models of Baugh et al. (1998). An angular 
fluctuation in the magnitude of the distortion would also 
be expected in this case, reflecting the discrete nature of 
the sources contributing to the effect. 

As seen above, reionization has several consequences; 
in the rest of this review we focus on the connection to 
the CMB, and the effects of inhomogeneities. 

3. EfTect on the CMB 

In this section, we give an in-depth treatment of the 
effects on the CMB photons of Thomson-scattering off of 
the free electrons produced by reionization. After defin- 
ing the optical depth and giving its dependence on the 
redshift of reionization and cosmological parameters, we 
discuss the three separate effects: damping, Doppler and 
polarization generation. Then we derive the evolution 
equations from the Boltzmann equation, argue that the 
Doppler effect is the most important inhomogeneous ef- 
fect, and calculate the resulting power spectrum for sev- 
eral simple models of inhomogeneous reionization (IHR) . 

The probability of a photon scattering in the time in- 
terval from some initial time U to the present, to is given 



by 1 — e '^'^*°^ where 



T(io) = / (JTne{t)dt 



(1) 



is the optical depth to Thomson scattering and ctt is the 
Thomson cross-section. If one assumes a step-function 
transition from a neutral to an ionized IGM at a redshift 
of Zion, the mean optical depth is given by (e.g., Grifhths 
et al., 1998) 



T(Zion) = T*/no 

where 



(l - Oo + fio (1 + Zion)^) 



1/2 



1 



Hoftbar 
47rG?Ti„ 



X (1 - Y/2) ~ omsnth, 



(2) 



(3) 



t{zu 



0.037/Vf^ 



l + Zj. 
11 



3/2 



0.02 



0.65 



(4) 
Thus we see that the fraction of CMB photons that have 
been scattered may be very small, but that this frac- 
tion grows fairly rapidly with increasing redshift. The 
increase in optical depth for fixed Zion with increasing 
1 — rig is due to the fact that the proper time from the 
present to Zjon increases as the cosmological constant in- 
creases. 

3.1. Simple Arguments 

Here we give simple arguments about the nature of 
the damping, Doppler and polarization effects. Before 
doing so, it is useful to define the two-point correlation 
function, C{9) which is perhaps the most important sta- 
tistical property of the CMB fluctuations: 

C{e) = (A(7i,x,?7o)A(72,x,ryo));cos6' = 71 -72 (5) 



where (...) denotes ensemble average. The same informa- 
tion is contained in its Legendre transform, the angular 
power spectrum, Ci: 



c{e) = J2 



2^ + 1 

An 



CiPeicosd). 



(6) 



fib and f2o are the density of baryonic matter and of all 
matter respectively, in units of the critical density today, 
Hq is the Hubble constant. 

Ho = 100/ikmsec~^Mpc~^, G is Newton's constant, irip 
is the proton mass and Y is the mass fraction of baryons 
in helium. The final equality assumes Y = 0.24. To get 
equation 0, we assumed that rig ex (1 -t- z)^, and adopted T -j- AT 

the line element, dz/dt, of a spatially flat universe with 
only non-relativistic matter and a cosmological constant. 
Note that for Zion = 5.64 (the observational lower 
bound), fifc/i^ = 0.02 and h = 0.65, we get r = 0.016 
for fio = 1 and r = 0.049 for ilp = 0.3. For Zion only 
slightly larger than this limiting value, it is a good ap- 
proximation to rewrite equation as 



Roughly speaking, Ci is the power in modes with wave- 
length tt/Z so that I = 180 corresponds to about a degree. 
If the fluctuations are Gaussian-distributed and statisti- 
cally isotropic, all other statistics can be derived from 
Ci. 

3.1.1. Damping 

For an initially homogeneous and isotropic photon 
field, there is just as much probability to scatter into 
a line of sight as there is to scatter out of it. Therefore, 
the optical depth to scattering produces no net effect. If, 
however, there is initial anisotropy, then this anisotropy 
is damped. Consider a line of sight along which the tem- 
perature differs by AT from the mean T in the absence 
of damping. If damping is present, then the temperature 
is changed to 



(T + AT) 
f + ATe~ 



(T + AT) (1 



T(l 



Rl 



1 - exp (-2t) 



a + cix + C2X^ + c^x^ + C4X^ 



+ exp(-2T), (8) 



(7) 



This equation expresses the fact that the final tempera- 
ture is given by the initial temperature, reduced by the 
photons that have been lost, and increased by the pho- 
tons that have been scattered in from other lines of sight. 
Since the photons that have been scattered in come from 
many different lines of sight, their average temperature 
is very nearly the global average — an assumption made 
in equation ^ The net result is that AT -^ ATe""^ and 
therefore C(0) -^ C(6')e-2^ or Ci -^ Cie~^^ . This simple 
calculation would imply that the damping is indepen- 
dent of scale, which is not exactly true. Our assumption 
that the average T along lines of sight that scattered in 
equals the global average holds only when the distance 
between the scattering surfaces is much larger than the 
length scale of interest. At very large angular scales, the 
damping does not occur, as one would expect from simple 
causality considerations. 

More precisely, the damping factor is /-dependent and, 
defined via C/ — Rf Cf"'^^'^"^ , is given by the fitting for- 
mula of Hu & White 1997: 



with X = l/{lr + 1) and a ^ -0.267, C2 = 0.581, C3 = 
—0.172 and C4 — —0.0312. The characteristic angular 



scale Ir is roughly the angular scale subtended by the 
horizon at the new last-scattering surface and is given 
approximately by (Griffiths et al. 1998, Hu & White 
1997): 



ir v-^ ^" ^ionj 



lony (1 



0.841nf7o)- 1. 



(9) 



If not for the ^-dependence, the damping effect would 
be completely degenerate with the amplitude of the pri- 
mary power spectrum; i.e., their effects would be indis- 
tinguishable. The similar response to amplitude and t 
over a large range of / makes them approximately degen- 
erate and is the reason why the optical depth can only 
be determined to about 10% (Zaldarriaga et al. 1997) or 
possibly worse (Eisenstein et al. 1998) based on temper- 
ature anisotropy alone. 

3.1.2. Polarization 

The Thomson scattering differential cross-section is 
polarization-dependent, and therefore the scattered ra- 
diation may be polarized even if the incident radiation 
is not. However, by symmetry considerations alone one 
can see that initially isotropic radiation will not become 
polarized^. In fact, it is easy to show that the particular 
anisotropy required for the creation of polarization is a 
quadrupole moment. This is essentially due to the cos 20 
dependence of the differential cross-section. 

The Fourier analogue of the angle- distance relation 
can be written as ^ ~ {rj — rj<f)k where k is the comoving 
wave-number that projects from conformal time r|.^, into 
multipole moment (. at conformal time rj. Conformal time 
is related to proper time by dr] = dt/ a where a is the scale 
factor of the expansion. Applying this relationship twice, 
and using the fact that polarization is generated by the 
quadrupole moment, ^ = 2, we find 

^p ^ {m - man) I ) - 2 i\fz~TT - i) . 

\Vion-VLSsJ 

. ^^^^ 

For the last equality, we assumed a matter-dominated 
Universe and that Tyion > > ilhss ■ Thus we expect a peak 
in the polarization power spectrum near ip ~ "^z^^^ . 

Determining the location of this peak allows one to 
determine the epoch of reionization. The amplitude of 
the peak is proportional to the optical depth, and thus 
also helps in determining Zion- Unfortunately, the signal 
is very weak and present only at large angular scales, 
e.g., €p ~ 5 for Zion = 10. Since there are only 21 + 1 
independent modes from which to determine Ci, sample 
variance is worst at low I. Also, the large-scale features of 



polarization maps reconstructed from time-ordered data 
are likely to be those most sensitive to systematic er- 
rors. Finally, polarized galactic foregrounds (dust and 
synchrotron) are expected to be most troublesome at 
small I (Bouchet et al. 1998, Knox 1998). However, a 
benefit of the low-£ location of this signature is that if it 
can be measured, its interpretation will not be compli- 
cated by the inhomogeneity of the reionization, which is 
a much smaller-scale phenomenon. 

If the systematic errors are negligible, the reionization 
feature in the polarization power spectrum can be used 
to discern very small optical depths. Forecasts of param- 
eter determination by Eisenstein et al. (1998) using both 
the temperature anisotropy and polarization data are for 
one sigma errors on r of 0.022 for MAP and 0.004 for 
Planck. Also see Zaldarriaga et al. (1997) for similar 
results. Note that these sensitivities, especially in the 
case of Planck, are high enough to detect the effects of 
reionization at any epoch given that we know it happened 
prior to z = 5.64. 

3.1.3. Doppler Effect 

The contribution to AT/T at location xq, in direction 
7, at conformal time 770 (today) is given by 

AT /"'" 

^p(xo,7, ?7o) = CTT / drjUe (x) J ■ Ve{x)a (11) 

where x = xo+7 (770 — rj). As one might expect, it is pro- 
portional to the line-of-sight integral of the parallel com- 
ponent of electron velocity, Vg times the number density 
of scatterers. Tig. 

As Sunyaev (1978) and Kaiser (1984) pointed out, in 
the homogeneous case this contribution to ^- is sup- 
pressed by cancellations due to oscillations in Vg. One 
way to think of these cancellations is that photons get 
nearly opposite Doppler shifts on different sides of a den- 
sity peak, which is a consequence of potential flows gen- 
erated by gravitational instability. The cancellations are 
less complete at large angular scales {I < 100) and there 
is a measurable effect which is exactly taken into account 
in Boltzmann codes such as CMBFAST^. If speed is crit- 
ical, the fitting formulae of Griffiths et al. (1998) can be 
employed for an approximate treatment. 

The cancellations at small scales can be greatly re- 
duced by the modulation of the number density of free 
electrons which we can write as Ue — Xefip where Xe is 
the ionization fraction and np is the number density of 
all protons, including those in H and He. This modula- 
tion occurs even in the case of homogeneous reionization 



'^Unless the electron spins are aligned, perhaps by a magnetic field. 



http://www.sns.ias.edu'^matiasz/CMBFAST/cmbfast.htm' 



(spatially constant Xg) due to spatial variations in Up. 
Because Ve is zero to zeroth order, and Up is uniform to 
zeroth order, this is a second-order effect. That this sec- 
ond order effect could dominate the first-order one (which 
is small due to cancellations) was pointed out by Ostriker 
and Vishniac and is known as the Ostriker- Vishniac ef- 
fect (Ostriker & Vishniac 1986, Vishniac 1987). Sub- 
sequent systematic study of all 2nd order contributions 
have shown that it is the dominant second-order contri- 
bution as well (Hu et al. 1994, Dodelson & Jubas 1995). 
Another source of modulation of the number density 
of free electrons is spatial variations in Xg, i.e., IHR. As a 
photon streams towards us from the last-scattering sur- 
face, it may pick up a Doppler "up kick" on one side of an 
overdensity, but avoid the canceling "down kick" on the 
other because the IGM there is still neutral. How effec- 
tively the cancellation is avoided depends on the match- 
ing between the typical sizes of the ionized domains and 
the correlation length of the velocity field, as we will dis- 
cuss below. 

3.2. Evolution Equations 

We now derive equation O, as well as a more general 
expression which includes the contribution from damping 
as well. More complete treatments of the evolution of a 
photon distribution function in an expanding, inhomoge- 
neous Universe can be found in, e.g., Ma & Bertschinger 
(1995). 

The Boltzmann equation governing the evolution of 
the photon phase-space distribution function, /, is sim- 
ply -Jf = C, where C takes into account the effect of 
collisions — in this case with electrons via Thomson scat- 
tering. The Thomson scattering cross-section is indepen- 
dent of frequency and therefore does not induce spectral 
distortions; the effect on the distribution function can 
be fully described by the change in brightness. Thus 
we use the brightness perturbation variable defined by 
A = {dfo/dq)^^fi where q is the comoving photon mo- 
mentum, and the phase-space distribution function has 
been partitioned into a homogeneous and inhomogeneous 
part, / = fo+fi- For a Planckian / with thermal fluctua- 
tions, AT, about a mean temperature of T, A = 4AT/T. 

Expanding the total derivative on /, keeping terms 
to first order in the perturbation variables and Fourier 
transforming the result, the Boltzmann equation be- 
comes: 



A(fe, 7, 7]) + ik^A{k, 7, 77) + -h{k, rf) 



drxe'^'^Tiyi, ?7)S'(x, 7, ??) 



(12) 



(3/i33(fc,ry)-/i(fe,ry))P2(Ai) 



where A is the Fourier transform of the brightness per- 
turbation, 7 is the direction of the photon momentum, 
fi = k ■ J, the dot indicates differentiation with respect 
to conformal time and ^,33 and h are the 3,3 component 
and trace of the synchronous gauge metric perturbation, 
respectively, in the coordinate system defined by x^ — 7. 
On the right-hand side of equation (^2|) we have explic- 
itly separated out the dependence of the collision term 
on the differential optical depth, f = aaq^ne (x, rj) . The 
quantity S* is a function of the brightness perturbation 
and the electron velocity. To first order in the pertur- 
bation variables its Fourier transform is (neglecting the 
polarization dependence of the Thomson cross-section) 

5(k, 7, 77) = -A + Ao + 47 • ve - iA2P2(M) (13) 

where A = ^^(2? + l)AePe{n). 

If we were to allow both f and S to depend on x, the 
Fourier modes would no longer evolve independently. To 
avoid this complication, we use an expansion in small f. 
Let A(") be the solution for f = and let A^^) = A-A(°). 
Then to first order in f: 

A^'\k)+7fcAiA(i)(k)= f £xe'^-''T{ji)S^''H^,l)- (14) 

which has the solution in real space: 

AW(7,x,,7o)= / d77'f(x')5('')(x',7) (15) 



where x' = x + 7(?yo ^ v')- This solution can be derived 
with a Green function approach or simply verified by 
substitution. It can be taken to higher order, if desired, 
by getting A*^^^ from S^^\ etc. Note that for the Doppler 
effect equation |l^ is exact because Ve (unlike, e.g.. A) 
receives no corrections due to the optical depth. 

We can now argue that the Doppler effect, due to the 
Vg term in equation WB, is the most important inhomo- 
geneous effect. The dominance of this effect is due to 
the fact that by the time of rcionization, velocities have 
grown substantially as they react to the gravitational 
field. On the relevant length scales, the rms peculiar 
velocities are about 10"'^. In contrast, the damping of 
the anisotropy is due to the —A term in S, which has 
not grown from its primordial value of 10~^. Similarly, 
polarization is sourced by the quadrupole moment which 
is also still at the 10~^ level at the epoch of reionization. 
Therefore we expect these contributions to be down by 
two orders of magnitude in amplitude, or four in power. 



If we break up the correlation function into compo- 
nents according to the expansion in optical depth: 

C{0) = C(") (9) + 2C("i) {9) + C(i) (9) (16) 

where C^^\ C^^^ and C^^^'> are the correlations between 
the two zeroth order terms, between the two first order 
corrections and between the first order and zeroth order 
terms, respectively. 

For the Doppler effect, the (exact) correction to the 
homogeneous case comes entirely from the C'^^-' term 
which is equal to 



C'^^\0) 



no 



drji 



Vo 



d?72((T7i •Ve)i(7-72 •Ve)2) (17) 



where the numeric subscript on the quantities in paren- 
theses means evaluate at, e.g., (x -I- 71(770 — 771), 771). 

From equation U^ we see that the relevant quantity 
for the CMB 2-point autocorrelation function, is the two- 
point function of the product, f7- Vg, integrated over the 
lines of sight. As pointed out by Knox et al. (1998) it 
is this latter 2-point function which must be calculated 
from either simulations or analytical models of the reion- 
ization process. 

3.3. Results from Simple Models 

A very simple toy model was used by Gruzinov and Hu 
(1998, hereafter GH) and by Knox et al. (1998, hereafter 
KSD) for investigating the effect of IHR on the power 
spectrum. In this model, independent sources turn on 
randomly and instantaneously ionize a sphere with co- 
moving radius R, which then remains ionized. The rate 
of source creation is such that the mean ionization frac- 
tion is zero prior to Zjon + Sz and then rises linearly with 
redshift to unity at Zion- GH derived the following ap- 
proximate expression for the resulting power spectrum: 



2tt 



R 



770 - '7ion 



A = 



36 



i<v'>R/ino6z{l + z,,^Y'' (18) 
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where {v^) is the variance of peculiar velocities and tq = 
cFTilono and hq is the number density of free electrons 
today. The qualitative shape of this power spectrum is 
simple to understand. On large scales it is the shape of 
a white noise power spectrum (C; = constant), due to 
the lack of correlations between the patches. At small 
scales there is an exponential cutoff corresponding to the 



angular extent of the patches themselves, since they have 
no internal structure. 

The exact calculations for this model by KSD result 
in a very similar shape but with the power reduced by 
a factor of two and the location of the peak shifted to 
larger I by 30%. Thus we can write down very simple, 
and accurate, equations for the location of the maximum, 
which are those of GH with the prefactors modified to fit 
the resuhs of KSD: 
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T^ <v^ > R/rjoSz (1 + Zionf^^ , and 



1 - (1 + z,,^r^/^ 



(19) 



Note that for sufficiently large i?, velocities are not 
coherent within a patch and equations |lg and |l9| break 
down; cancellations once again become important. How- 
ever, this only happens for i? > 10 Mpc since the velocity 
correlation function only becomes negative for separa- 
tions greater than about 30/i~^ Mpc (KSD). As long a 
it! is less than lO's of Mpc, the power scales linearly with 
R. 

Aghanim et al. (1996, hereafter "A96") used a slightly 
more complicated model in which there was a distri- 
bution of patch sizes due to an assumed distribution 
of mini-quasar luminosities. The larger patches domi- 
nate at low I and the cutoff at higher I is much softer 
than in the single patch case, due to the existence of the 
smaller patches. For an order-of-magnitude calculation, 
their largest patches have i? ~ 10 Mpc and reionization 
proceeds from Zion + <5z = 10.6 to zion = 5.6. For this sim- 
plification of the A96 model we find the curve in figure @ 
which has (^) "^l^ lO^^^ ^^d /^ax '^ 1300. 

V / max 

This amplitude is just large enough to possibly be a 
source of systematic error in the determination of cosmo- 
logical parameters from MAP data, which has sensitivity 
out to I ~ 1000. To be more precise, MAP can deter- 
mine the power in a broad band centered on I = 800 
with width &l = 400 to about 0.5 %. Thus, given the pri- 
mary signal assumed in the figure, we may expect signals 
on the order of 5 x 10~^^ to be significant contaminants. 
This IHR signal from a 10 Mpc patch would be an even 
more noticeable effect in the Planck data since Planck 
has sensitivity out to about I ~ 2000 and can determine 
the power in a broad band centered on / = 1500 with 
width 51 = 1000 to about 0.2 %. 

The quasar luminosity function (QLF) is relatively 
well measured in the optical below redshift 2: < 5, with a 
typical average optical luminosity of 10''^ erg s~^. In 
order to estimate the expected typical HII patch size 
around the observed quasars, one needs to know their 
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Fig. 3. — The thick sohd hne is the angular power spec- 
trum of the primary CMB anisotropy in the standard 
cold dark matter (SCDM) model, normalized so that 
the fractional rms fluctuations of the mass in 8/i^^ Mpc 
spheres is a% — 1.2. The other curves are contributions 
to the power spectrum expected from different models 
of IHR. The light solid curves are for the uncorrelated 
model (see equation |l8| ) with Zjon = 5.6 and 5z = 5 
with comoving patch radii of _R = 10 Mpc and i? = 1 
Mpc. The dashed curve is for the correlated model of 
KSD which achieves full reionization at Zjon = 26. The 
underlying model is spatially flat CDM with fif, = 0.1 
and h = 0.5. The expected 0-V effect for this same 
model with the same reionization redshift is also shown, 
calculated using the formalism of Jaffe & Kamionkowski 
(1998). 



lifetime (or "duty-cycle"). Taking the lifetime of quasars 
to be universal, and near the Eddington time (~ 5 x 10^ 
yr), and further ignoring recombinations, one obtains a 
comoving radius of ~10-20Mpc. As seen above, this 
would be large enough for its effects on the CMB to be 
detectable by Planck and maybe MAP. 

We emphasize however, that the typical quasar lumi- 
nosity at z = 10 can be much lower than at z ~ 3. 
In hierarchical structure formation models, the nonlinear 
mass-scale is ~ 3 orders of magnitude smaller at 2 = 10 
than at z ~ 3. With a linear scaling between quasar lu- 
minosity and halo mass, this would reduce the average 
patch-size by a factor of ~ 10. In addition, the duty- 
cycles of quasars are not well established. In theoretical 
models where the QLF is derived by associating quasars 
with collapsed dark halos (in the Press-Schechter formal- 
ism), a typical lifetime of '^ 5 x 10^ yr would generically 
lead to an overestimate of the number of observed quasars 
(e.g. Smallfe Blandford 1992). To avoid these overpre- 
dictions, one needs to assume that either (1) the quasar 
shines at a small fraction of the Eddington luminosity, 
or (2) the efficiency of black hole formation, expressed 
as Mbh/Mhaio, is 2 orders of magnitude lower than sug- 
gested by nearby galaxies (Magorrian et al. 1998), or 
(3) the black hole masses grow in an optically inactive 
phase (e.g. Haehnelt, Natarajan & Rees 1998). Alterna- 
tively, one can fit the observed QLF by assuming that the 
quasar lifetime is ~100 times shorter than the Eddington 
time (Haiman & Loeb 1997). In the latter case, the HII 
patch sizes would be reduced to i? ^ 1 Mpc (i.e. reion- 
ization would correspond to a larger number of smaller 
patches) . 

According to equation |l9|, patches smaller than a few 
Mpc result in power spectra with smaller maxima at val- 
ues of I that are too high, and well outside the range of 
sensitivity of MAP and Planck. Even for patches as large 
as 1 Mpc, one can see that, at least for reionization at 
the extreme low end of the redshift range, their effect is 
negligible (see figure ^. A detection of an IHR signal 
in the CMB with a peak could have important implica- 
tions for quasar evolution theories, e.g. by ruling out 
models that lead to short duty cycles and small patch 
sizes. A detection might also prove that reionization was 
caused by mini-quasars, rather than mini-galaxies, since 
for the mini-galaxy models, the patch sizes are expected 
to be only a few hundred kpc. However, in light of the 
small expected patch sizes, it seems unlikely that MAP 
or Planck reaches the required angular resolution for such 
a detection. 

The above calculations of the IHR power spectrum all 
assumed that the ionized regions are uncorrelated. How- 
ever, the pattern of density fluctuations will affect the 
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pattern of ionized regions and therefore correlations in 
the density field should give rise to correlations in the 
ionization fraction field. Determining what these corre- 
lations should be is complicated by the fact that the den- 
sity field affects the ionization field in two different ways: 
by the fact that the sources are most likely located in 
high density regions and the fact that the recombination 
rates are highest in the densest regions. 

To take the former effect into account, KDS gave the 
ionization field the same correlation structure as that 
of the "high-peaks" of the density field. These high- 
peaks are regions where the fractional density pertur- 
bation, 5 = —, smoothed on the appropriate scale, 
exceeds the critical value for spherical collapse. The 
smoothing scale was taken to be that corresponding to 
1O®M0[(1 -I- 2;)/ll]~'^/^, since objects below this size will 
not be able to cool sufficiently to fragment and form stars 
(Haiman, Rees & Loeb 1997). One further assumption 
was that the sources ionized a region E times larger than 
that from which they collapsed where the efficiency factor 
is time-dependent and given by (Haiman & Loeb 1997, 
Haiman & Loeb 1998a) E = 7 x 10^{r]/Tjof. For this 
rather high choice for the efficiency, reionization is com- 
plete at Zion = 26. 

The resulting spectrum is shown as the dashed line in 
figure g. The correlations between the patches — which 
are individually only hundreds of kiloparsecs in radius — 
drastically alter the low I shape. Even given the high 
redshift of reionization in this scenario, the corresponding 
uncorrelated model of equation |l8| would predict orders 
of magnitude less power at I ^ 10^ to 10*. 

If this signal were in the data, but not in ones model of 
the data, it would lead to a bias in the determination of 
cosmological parameters. KSD estimated the bias that 
would occur in a nine parameter fit, assuming the pri- 
mary signal was that of standard cold dark matter. For 
Planck, the systematic in the estimate of the baryon den- 
sity was equal to the statistical error. For all the other 
parameters the systematic error was less than the statis- 
tical error. 

3.4. Discussion 

If the mini-quasar scenario is correct, with bright, 
long-lived quasars, then the Doppler effect will lead to 
a significant contamination of the Planck data and pos- 
sibly a marginally important one for MAP data. For 
the more likely mini-galaxy scenario, if one ignores the 
correlations in the ionized patches, the contamination is 
completely irrelevant for both MAP and Planck. How- 
ever, correlations may indeed be important, and in one 
attempt to take them into account, we have seen that 



for an extreme mini-galaxy scenario, the contamination 
is marginally significant for Planck. 

Work by Oh (1999) supports the idea that the ion- 
ized patches should be correlated. Oh has modeled the 
ionizing sources as characterized by a mean separation 
(given by Press-Schechtcr theory), an isotropic attenua- 
tion length which specifies how far the ionizing radiation 
propagates, and a power-law correlation function (appro- 
priate for these highly-biased objects). Placing sources in 
a box and working out the resulting fiux densities he can 
then calculate the correlation function of the flux density. 
It is indeed correlated, with significant (0.2) correlation 
at a comoving separations of 10 Mpc. Although a clump- 
ing factor for the IGM from Gnedin & Ostriker (1997) is 
used to estimate the attenuation length, all propagation 
from the sources is taken to be isotropic, and hence this 
work does not yet include the influence of shielding and 
shadowing. 

As one can see in Fig. H, the contribution from IHR 
may be sub-dominant to that from the Ostriker- Vishniac 
(0-V) effect, as it is for all three models of IHR shown. 
The 0-V contribution here is calculated using the formal- 
ism of Jaffe & Kamionkowski (1998, hereafter "JK98"). 
One must keep in mind that as a second-order effect, it 
is highly sensitive to the normalization. The curve here 
is COBE- normalized (as — 1-2); a cluster normalization 
of (Tg = 0.6 would reduce the Ci by a factor of 16. It 
is also calculated for the very early reionization redshift 
(and resulting high optical depth) of zion — 26. However, 
the amplitude is not very sensitive to r because much 
of the effect actually comes from more recent times — the 
dropping number density of electrons must compete with 
the growth in the density contrasts and peculiar veloci- 
ties(Hu & White 1995, JK98). Of course, the separation 
into IHR and 0-V is artificial, and efforts to understand 
the contribution from reionization are surely going to 
benefit from treating them simultaneously. Indeed, the 
2-point function of f-y • Vg, measured from a simulation, 
would include both effects. 

Gravity is the only other late-time, frequency-independent, 
infiuence on the CMB photons, causing lensing and the 
Rees-Sciama effect (Rees & Sciama 1968). The angular- 
power spectrum from the Rees-Sciama effect has been 
calculated by Seljak (1995) and is substantially sub- 
dominant to either the primary spectrum or the 0-V ef- 
fect at all angular scales. Gravitational lensing results in 
a smearing of the primary CMB angular power spectrum 
and has also been calculated by Seljak (1996). Uncer- 
tainties in this calculation may indeed be important for 
any attempt to recover the contribution from the IHR 
and Ostriker- Vishniac effects. 
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4. Conclusions 

The principal unanswered questions about reioniza- 
tion are: what type of sources caused it, and at what red- 
shift did it occur? These questions have been addressed 
both theoretically and observationally. From a theoreti- 
cal point of view, the two leading candidate sources are 
an early generation of stars ( "mini-galaxies" ) , or massive 
black holes in small halos ("mini-quasars"). It is possi- 
ble to estimate the efficiencies with which objects form in 
the earliest collapsed halos: from the metals in the Lya 
forest, and from the evolution of the quasar luminosity 
function at z < 5, respectively. However, the uncertain- 
ties in these efficiencies are still too large to allow definite 
predictions. The expected reionization redshift, depend- 
ing on the type of source, cosmology, power spectrum, 
and a combination of the efficiency factors, is between 
7 ^ z ^ 20 in homogeneous models. A further compli- 
cation is the inhomogeneous nature of reionization: the 
sources of ionizing radiation are likely clustered, and em- 
bedded in complex dense regions, such as the filaments 
and sheets seen in 3D simulations. At the same time, 
the gas to be ionized likely has significant density fluctu- 
ations. 

Further theoretical progress on the problem will likely 
come from 3D simulations. Such simulations must be 
able to follow the propagation of a non-spherical ioniza- 
tion front into a medium that has large density fluctu- 
ations, as well as opaque clumps of absorbing material. 
In a homogeneous medium populated by randomly dis- 
tributed sources, the background radiation would have 
negligible fiuctuations (e.g. Zuo 1992). However, be- 
cause the density inhomogeneities absorb the UV flux of 
the ionizing sources, and re-emit diffuse ionizing photons 
as they recombine, the ionizing background will likely be- 
come significantly inhomogeneous and anisotropic (Nor- 
man et al. 1998). The first numerical algorithms to deal 
in full detail with these problems have been proposed by 
Abel, Norman & Madau (1998) and Razoumov & Scott 
(1998) (see also Gnedin 1998). 

The reionization redshift can also be constrained ob- 
servationally, with a resulting uncertainty that is com- 

rablp to the thenretira] rangp- P> < y < 40 C)n both 



para 

ends of this redshift range, nbservationa.1 progress is likely 
in the near feature. If new Lya emitters are discov- 
ered at higher redshifts, this would improve the lower 
bound. Conversely, as the CMB data are collected, the 
constraint on the electron scattering optical depth will 
improve, tightening the upper bound. 

We have seen how the optical depth to Thomson- 
scattering, and hence the epoch of reionization, can be 
constrained by determining the amount of damping in the 



CMB temperature anisotropy and by detection of the po- 
larization contribution at very large angular scales. The 
spatial inhomogeneity of the reionization process gives us 
a further opportunity to probe in more detail the epoch 
of reionization. For this probe to become reality, fur- 
ther development of theoretical predictions is required, 
as well as more sensitive measurements at arc minute 
and smaller angular scales. 

Fortunately, the inhomogeneity of reionization is un- 
likely to spoil our ability to interpret the primary CMB 
anisotropy, to be measured with exquisite precision over 
the next decade, although this may not be the case if 
reionization did indeed occur via mini-quasars. For the 
mini-galaxy models, even when reionization occurs at the 
very high end of the allowed redshift range, very little sig- 
nal is produced at the relevant angular scales. Further 
progress in these predictions is likely to be motivated by 
the desire to understand the spectrum at I > 3000 as a 
means of probing the end of the dark ages, rather than 
at I < 3000, as a possible contaminant of the primary 
signal. 

We are grateful to our collaborators S. Dodelson, A. 
Loeb, M. Rees, and R. Scoccimarro for contributing to 
our understanding of this subject, to A. Jaffe for supply- 
ing the "0-V" effect curve in Fig. |^, and to M. Norman 
for supplying Fig. 2. ZH was supported at Fermilab by 
the DOE and the NASA grant NAG 5-7092. 
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